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The effect of steady grazing flow on the liner is included in the 
model. Sample calculations of liner acoustic resistance with spectral 
excitation both with and without grazing flow are presented. 
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SYMBOL LIST 

th 2 

amplitude of K frequency component of P , lbf/ft 

resonator cavity depth, ft 

amplitude of K th frequency component of v , ft/sec 
speed of sound, ft/sec 

coefficient of Fourier cosine series fit to v , ft/sec 

o 7 

orifice diameter, ft 
orifice effective length, ft 

orifice length (perforated plate thickness), ft 
grazing flow Mach number 

2 

excitation pressure at orifice, lbf/ft 

th ■ 2 

K frequency component of P q (peak amplitude) , lbf/ft 

grazing flow effect on specific acoustic resistance of an orifice 
(eq. (11)) 

2 

total orifice resistance, lbm/ft sec 

nonlinear orifice resistance, lbm/ft sec 

sound pressure level (ref. 2 xlO~^ microbars) 

coefficient of Fourier sine series fit to v , ft/ sec 

o 

time, sec 

initial and final time over which Fourier series fit to v is 
made (see eqs. (19) and (20)), sec ° 

orifice fluid particle velocity, ft/sec 

root-mean-square orifice particle velocity, ft/sec 

l 

steady flow grazing velocity past liner, ft/sec 
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V 
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maximum fluid displacement in orifice, ft 

fluid displacement in orifice, ft 
th 

phase angle for K frequency component of P , rad 
orifice end correction, ft 

specific acoustic resistance of an array of resonators 

specific acoustic resistance in orifice (0 = g0) 

frequency dependent part of 0 QL (eqs. (15) and (16)) 

specific acoustic resistance of the K frequency component of P 
in the orifice ° 

linear specific acoustic resistance in orifice 

nonlinear specific acoustic resistance in orifice 

frequency independent part of 0 (eq. (14)) 

UJL 

specific acoustic resistance of an array of resonators at the 

fundamental frequency when the harmonic is at much lower amplitude 

dynamic gas viscosity, lbm/ft-sec 

kinematic gas viscosity, ft /sec 

3 

gas density, lbm/ft 


o resonator array open area ratio (orifice area/wall area) 

4> k phase angle of K frequency component of v , rad 

X specific acoustic reactance of an array of resonators 

Xq specific acoustic reactance in orifice (x - ox) 

Xqk specific acoustic reactance of the frequency component of P 

in the orifice 

a) angular frequency, rad/sec 


RESONATOR IMPEDANCE MODEL 

The geometry of the model and the important dimensions are shown in 
figure 1, The model which will be presented approximates the actual 
fluid motion with simplified one-dimensional flow equations- 
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The fluid motion between stations 0 and 1 is considered to be only 
mass reactance (inductive) and resistance controlled. Some of the dissi- 
pation actually occurs outside of the orifice in the jets formed by the 
flow (ref. 3), but the resistance is considered in the equations of the 
orifice since nonlinear resistance can be related to orifice velocity. 

The mass continuity equation describes the conditions between sta- 
tions 1 and 2. Where a relationship between pressure and density is re- 
quired, the isentropic relation is used. For the calculations presented 
in this paper the back cavity (2 to 3) was divided into four equal seg- 
ments with both mass reactance (inductance) and stiffness reactance 
(capacitance) considered. In this way the acoustic transmission line 
properties of the back cavity could be approximated without actually cal- 
culating a finite difference solution using incremental cavity lengths. 

In the equations which follow, however, lumped stiffness in the cavity is 
assumed to keep the equations simple for purposes of illustration. 


When the system of equations 
dX 


d 2 X 


PL 


+ R 


dt 
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are combined the results are: 


P 
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and 

dX 

v « — — 

o dt 

The effective length (L) of the orifice is given by: 
L = £ + e 

where e is the orifice end correction (ref. 11): 
0.85 d (1 - 0.7 Vo) 

e 3 

1 + 305M 


( 1 ) 

( 2 ) 

( 3 ) 

(4) 


The coefficient of the second derivative is associated with the induc- 
tance and that of the first derivative with the resistance of the system. 
The capacitance is associated with the coefficient of the displacement. 

The subscript on X and P denotes that they are the displacement and 
pressure in the orifice. 

When solutions of equation (1) are considered it should be noted, as 
previously stated, that the more complicated system approximating wave 
travel in the back cavity will be used. 

If the resistance R Q in equation (1) were a constant, the pressure 
could be given as a series of sinusoidal functions and a closed form solu- 
tion would result. However, it is of interest here to consider a nonlinear 
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resistance term which will be considered to be a function of the orifice 
velocity (v Q ) and displacement (X Q ) as well as the grazing flow velocity 
( V oo) • The solution must thus be obtained by numerical integration. 

The definition of the nonlinear resistance will now be developed. 


Effect of Orifice and Grazing Plow Velocities on Nonlinear Resistance 

The nonlinear resistance of an orifice has been shown to depend upon 
the magnitude of the orifice velocity (refs. 4, 5, and 6). This resist- 
ance can be written as: 

R 0N = p l v ol (5) 

A proportionality constant is often used to account for the effects of a 
discharge coefficient; however, this will not be used here. 

It is convenient to use the normalized or specific resistance defined 

as : 


9 0N pC C (6) 

which is just the magnitude of the Mach number in the orifice. 

The influence of a grazing flow velocity on steady flow orifice re- 
sistance can be seen in figure 2 in which the steady flow resistances of 
several orifices are presented both with and without grazing flow. This 
data is presented through permission of Pratt and Whitney Aircraft and 
the Boeing Company who jointly obtained this data. The steady flow re- 
sistance data will be used to develop an acoustic resistance model by 
equating the instantaneous acoustic resistance to the steady flow resist- 
ance. 


It can be seen in figure 2 that when v m = 0 the data can be approx- 
imated by the solid lines which are given by equation (6) . This is true 
for both positive and negative orifice velocity. However, when a grazing 
flow velocity is present the results change substantially. 

Three regimes can be identified which depend upon the sign and mag- 
nitude of the orifice velocity (v Q ) . When v Q >_ 0, e 0N appears to be 
increased by a constant amount which is proportional to v^,. This remains 
true for all values of positive v Q shown. For negative v Q , however, 
there are two regimes. There is a transition region in which 6 QN is 
approximately constant with v Q and proportional to v^. At sufficiently 
large negative values of v Q , 0 q N is again approximated by the solid 
line given by equation (6) . 

The equations used to describe these regions can be expressed as 
follows: 
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If v > 0 
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6 ON 
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where 


Q = 



if v < 0 
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6 0N 


- Q 




6 ON 





(7) 


( 8 ) 


(9) 

( 10 ) 


As can be seen from the data In figure 2 the intercept of the dashed 
lines with the ordinate at v Q = 0 should be considered to be a function 
of open area ratio (a)- The definition of Q in equation (8) was thus 
further modified to be: 


Q = 0.3 ^ ^ 3 ' 15 * 10 + 1.77 a + 0.3^ (11) 

The relationships of equations (7), (9), and (10) were assumed to 
hold instantaneously and modified forms to account for orifice displace- 
ment were used to define the nonlinear resistance needed in equation (1). 

The modified forms will now be presented. 


The Effect of Displacement in the Orifice on Nonlinear Resistance 

Measurements of acoustic impedance of perforated plates with single 
frequency sinusoidal excitation show that the resistance given by equa- 
tion (6) is not attained at high frequencies. Although the fluid velocity 
Is sufficiently high to provide a nonlinear resistance, the fluid dis- 
placements are low at high frequency. 

The modification of the nonlinear resistance to account for fluid 
displacement might be justified by the following argument. The nonlinear 
dissipation is usually attributed to the energy lost in the jets formed 
by high amplitude oscillations (ref. 3). Unless the displacement of 
fluid within the orifice is sufficiently large the jet will not fully 
form and the amount of nonlinear dissipation (and thus resistance) will 
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be reduced e It was thus assumed that the fluid must be displaced a por- 
tion of the orifice diameter before a jet is fully formed. The axial 
extent of the turbulent mixing region can be related to the orifice diam- 
eter and thus provides the incentive to use the diameter as a normalizing 
factor . 


The modification of equation (6) required to agree with experimental 
results is given by the following equation. 



( 12 ) 


Equation (12) is applied instantaneously in the numerical integration. 

The quantity ^ is the maximum displacement obtained during a cycle. 
This is the displacement obtained when the orifice slug comes to rest and 
starts to move in the opposite direction. The exponential in equa- 
tion (12) is seen to act as an on-off function applied to the nonlinear 
resistance. When the displacement from maximum is much less than the 
orifice diameter (d) , then 0p N = 0. For - X Q ) » d, 0 ON * (vJ/C. 
The exponential in equation (12) was derived from one set of resistance 
versus frequency data. The justification for its use is provided by an 
adequate fit to a wide range of data. 

It should be noted that the orifice fluid displacement correction 
was not applied to the grazing flow term (Q) in equations (7) to (10) but 
only to the |v Q |/C terms. 


Linear Resistance 


For very small pressure amplitudes the acoustic resistance is deter- 
mined by the viscous losses in the orifice. This may be accounted for by: 


OL 


- 0 


OS 


+ 0 


of 


(13) 


The steady flow resistance 0Qg is included to provide a low frequency 
limit for © and is given by: 

0QS = 5.4x10^ y£ (14) 

where i is in feet. Equation (14) was derived from a correlation given 
in reference 7. 


For single frequency sinusoidal excitation the frequency dependent 


component of e is given by 


°of 




(1 + i/d) 


(15) 
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When multiple frequency excitation is used a derived form of equation (15) 
can be given by 



(16) 


Where the P are the several peak pressure amplitudes. Equation (16) 
is not a particularly satisfactory expression for the linear resistance, 
but it is all that was available. What is needed is an approximate solu- 
tion to the boundary layer equations of a flat plate with arbitrary motion 
which would replace the sinusoidal motion result of equation (16) . For- 
tunately, the linear resistance terms are usually insignificant when com- 
pared to the nonlinear resistance terms especially when grazing flow is 
considered. 


Equation (13) (using eq. (14) and (16)) is added to the nonlinear 
resistance to provide the total resistance for use in equation (1) . 


CALCULATION PROCEDURE 


The forcing function (pressure P Q ) in equation (1) can be an arbi- 
trary function of time. However, it will be useful for the acoustic im- 
pedance calculations to describe P by: 


P 

o 




sinCw^t + 


V 


(17) 


where 03 ^ and are chosen to simulate the desired acoustic power 

spectrum. Since phase information is usually not available, the phase 
angles (6^) are selected randomly. 

Equations (1) and (2) are then numerically integrated using the 
Runge-Kutta method (ref. 8). The initial conditions are that velocity 
and displacement are zero. The numerical solution is allowed to progress 
for sufficient time for the transient solution to decay. A Fourier 
series fit to the calculated velocity is then started and subsequently 
monitored at multiples of the period of the lowest frequency component. 
The orifice velocity can be expressed as: 


v o (t) 


( C K COS “k*" + S K S * n 


K 


(18) 
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with the same frequency components considered as in equation (17). The 
Fourier series fit over an arbitrary time internal with frequency compo- 
nents which are not necessarily multiples of a fundamental is more com- 
plicated than the standard series fit. The results are the same, how- 
ever, if a sufficiently long time interval is considered. This can be 
shown by the following expressions for the coefficients: 


C K ~ 


v (t)cos to t dt + F(ui.t,C . ,S .) 
, t ° K J 3 Y 

_1 

(t 2 - V 


( 19 ) 


and 


S K = 


v (t)sin to t dt + G(to . t ,C . ,S . ) 
/r ° K J J 3 

1 

(t 2 - t 1 ) 


( 20 ) 


The ratios of the integrals to (t 2 ~ t-^) will reach finite limits as 
( 1 2 “ The functions F and G, however, are purely oscillatory 

in nature and the limit of F or G divided by (t£ - t|) will approach 
zero. 


Once the Fourier coefficients C K and S K are determined the ve- 
locity can be expressed in a more convenient form as: 


V o (t > 



sin(oj t + 

is. 


♦k> 


( 21 ) 


which is an alternate form of equation (18) with phase and amplitude 
given by: 


<t> K = arc 



( 22 ) 


and 


B 


K 




( 23 ) 


Using the familiar equation for acoustic impedance as pressure divided by 
velocity there results 


0 




OK pCB. 


K 


cos(6 k - (f) K ) 


(24) 



and 
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X OK = pCB sin(6 K <t> K ) - < 25 ) 

K 

Equations (24) and (25) express the specific acoustic resistance and reac- 
tance for the frequency component in the resonator orifice. These can 
be converted to resonator array values for comparison with experimental 
data by dividing each by the open area ratio (a) of the array. 


COMPARISON WITH EXPERIMENTS 

This model was developed to provide acoustic impedance calculations 
when the exciting pressure contains a multitude of frequency components. 
Unfortunately no such data exists at this time for a direct comparison. 
There is a small amount of two frequency data available and a comparison 
with this data will be made. Single frequency data currently exists which 
will be used to check the analytical model for the effects of orifice 
fluid velocity and displacement and also grazing flow velocity. 


Single Frequency Without Grazing Flow 

Specific acoustic resistance data obtained by Garrison (ref. 9) using 
the two-pressure measurement method are shown in figure 3. The data are 
for four samples at several frequencies with a constant sound pressure 
level (160 dB) . The data represented by circles (£ = 0.05 in.) were used 
to derive the fluid displacement dependence expressed in equation (12). 

Two calculated resistance curves are shown for this condition at high 
frequencies with and without the displacement correction. The need for 
such a correction is apparent from figure 3. The comparisons with other 
orifice lengths (sheet thickness £) is also seen to be good. The maxi- 
mum errors occur with £ = 0.2 inch for which the data appears to be out 
of progression with the rest of the data at the lower frequencies. 

More data from Garrison (ref. 9) are shown in figure 4. A wide 
range of open area ratios (a = 0.015 to 0.13) and sound pressure levels 
(SPL = 130 to 168 dB) are included. The comparisons with the calculated 
resistances are seen to be good. Where differences occur there does not 
seem to be a consistent pattern with open area ratio or sound pressure 
level and these errors may well be within the accuracy of such measure- 
ments . 

In figure 5 the specific acoustic reactances are shown for four 
samples (ref. 9). Again the agreement with the calculated results is 
considered adequate. 
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Single Frequency with Grazing Flow 

Specific acoustic impedance measurements made by Garrison (ref. 9) 
are shown in figure 6. Two liner samples at two sound pressure levels 
were investigated using the two pressure method with grazing flow. The 
data at v^ =0 are questionable since the electro-pneumatic driver 
provided substantial air flow which had to be bled off before reaching 
the sample and circulating flows resulted (conversation with Garrison) . 

For both samples, at the lower sound pressure level (150 dB) and 
high velocities, the resistance data is increasing with v** indicating 
that the grazing flow effects are beginning to dominate the resistance- 
In this region the calculated results are in excellent agreement with 
the data. For lower values of v^ (especially for the a = 0,0272 
sample at 150 dB and the a = 0.0212 sample at 160 dB) the calculated 
resistance underpredicts the data which are virtually flat with changing 

v . 

00 


A possible explanation for the deviation between the data and the 
calculations at low grazing flow velocities can be made. Only in the 
case of a = 0.0272 (160 dB) did the orifice effective length, inferred 
from the data, resemble that of the model (eqs, (3) and (4)). For the 
other data the effective length decreased much more rapidly with v ro * 

In fact it dropped below the physical orifice length (£) which is not 
allowed by the present model. The reactance was thus dropped from a 
fairly high positive value toward the sample tuned point (x = 0) . The 
nonlinear resistance would then increase since jv Q | is roughly propor- 
tional to (6 Q + X Q ) ' . The current model is unable to account for this 
large change in orifice effective length. This is an apparent weakness 
in the model which will be resolved only with the testing of a more re- 
fined instantaneous effective orifice length assumption and perhaps the 
study of more data. 

In the regime of lower sound pressure levels and high v^ (of prac- 
tical interest in turbofan suppressors) the agreement between the data 
and the theory as shown in figure 6 is quite promising. 


Two Frequency Results 

The results of two frequency resistance measurements made by Garrison 
(ref c 1) for two liner samples are shown in figure 7. This data was taken 
as follows. The amplitude of the fundamental frequency (2000 or 2200 Hz) 
was held constant at 135 decibels. The amplitude of the harmonic was 
then varied from 110 to 150 decibels and the impedance at the fundamental 
was measured. The resistance ratio was then formed with 0 RE p being the 
resistance at the fundamental frequency for low harmonic amplitude 
(*110 dB). 

Calculated resistance values were obtained under the same conditions. 
Qualitative agreement is seen between the theory and the data in that the 
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second harmonic amplitude, at high amplitude, determines the resistance 
at the fundamental frequency. To produce agreement between theory and 
data, the second harmonic amplitude must be about 4 or 5 decibels higher 
in the calculation for one sample (a = 0.054) and about 10 dB higher for 
the other (a = 0.08). It is also disappointing that the trends (ampli- 
fication of resistance by the harmonic) between the two samples is oppo- 
site in the data to that in the theory. 

It is not obvious what changes in the theory would be required to 
get agreement with the two frequency data without destroying the excel- 
lent agreement that was obtained with the single frequency data. It is 
also not obvious that the theory should be altered on the basis of this 
small amount of data. 

A calculated result of Groeneweg (ref. 10) is also shown in figure 7. 
This uses overall sound pressure level for the two frequencies applied to 
a single frequency model as in reference 11. This method overpredicts 
the interaction effect since it essentially applies the entire harmonic 
pressure amplitude to the fundamental frequency impedance calculation. 

The agreement between theory and experiment is not at this time 
completely satisfactory. However, it is not clear which is at fault 
until additional data is thoroughly investigated. The interaction effect 
of one frequency upon another is clearly demonstrated to be of importance 
by both the theory and experimental data. 


SAMPLE CALCULATIONS WITH SPECTRAL EXCITATION 

The exciting pressure spectra used in the sample calculations are 
shown in figure 8(a). Six pressure terms (for use in eq. (1)) are used 
with a constant sound pressure level of 130 decibels. These are spaced 
at 500 Hz intervals and are intended to provide a crude simulation of 
white noise. The phases are chosen randomly. At 1875 Hz a pressure com- 
ponent is used with its amplitude varying from 130 to 150 decibels. This 
component is intended to simulate a discrete tone standing at or above the 
white noise. Hereafter these simulations will be referred to as the white 
noise and the discrete tone. 

In figure 8(b) some sample calculations of specific acoustic resist- 
ance (0) with zero grazing flow velocity (v ot = 0) are shown. The inter- 
action effect of the discrete tone with the background white noise is 
clearly demonstrated. When the discrete tone is at the same level as the 
white noise (130 dB) , the resistance spectrum is fairly flat with values 
around 8 = 0.25. When the discrete tone is increased to 140 dB the re- 
sistance at the white noise frequencies increases to about 0.5 and at 
150 dB to about 1.1. An interesting observation in figure 8(b) is that 
the resistances at the white noise frequencies are increased more than 
the resistance at the discrete tone frequency itself. This appears to 
be a general trend for the limited number of sample calculations com- 
pleted thus far. 
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In figure 9 the same pressure spectrum as shown in figure 8(a) is 
applied to the acoustic liner but this time with a grazing flow velocity 
v oo = 400 feet per second. Here the resistance is mainly grazing flow 
dominated and the amplitude of the discrete tone has much less effect. 
Only at the highest discrete tone amplitude (150 dB) are the resistances 
at the white noise frequencies increased and then only by a small amount. 
Again the resistances at the white noise frequencies are increased more 
than at the discrete tone. 

Small irregularities in 0 with frequency should not be considered 
as significant since the resistances are somewhat effected by the choice 
of pressure phase angles. This is demonstrated In figure 10 where three 
choices of phase angles were used both with and without grazing flow. 

The phase angles were still chosen at random, but three different random 
sets were used. The variations with phase angle are not large, but per- 
haps an average should be used over several sets of phase angle choices 
to better define the resistance. This will have to be done since phase 
information is usually lacking. 


CONCLUDING REMARKS 

Several investigators have used a single frequency sinusoidal exci- 
tation impedance model in an attempt to model an acoustic linear with 
spectral excitation (refs. 11 to 13). Two methods have been used to mod- 
ify these models to account for the multiple frequency excitation. 

The first method uses the overall sound pressure level of the noise 
spectrum (refs. 11 and 12). This results in applying the same pressure 
amplitude at each frequency. Since the reactance varies with frequency 
and helps determine the peak orifice particle velocity, the nonlinear re- 
sistance at each frequency will be different. This method results in a 
low resistance at both high and low frequencies with a maximum at the 
liner tuned frequency (X = 0) . 

The second method as used in reference 13 and discussed In refer- 
ence 11 uses a root-mean-square particle velocity (v^g) approach. A 
peak orifice particle velocity is calculated at several frequencies using 
the sound pressure level at each frequency. Then v^g is formed from 
the several velocities to give a single velocity to represent the spec- 
trum, This root-mean— square velocity is used in an equation similar to 
equation (5) (or an equation derived from expressions like eqs . (7) to 
to (10) if v^ were considered) to calculate resistance. Thus the re- 
sistance spectrum is uniform with frequency. 

From the preliminary calculations presented In this paper in fig- 
ures 8(b), 9, and 10 the resistance spectrum appears to be fairly uniform 
with the possible exception of the frequencies at which the high ampli- 
tude discrete tone appears. This uniformity suggests that the v i?mc 
method is at least qualitatively more correct than the overall sound 
pressure level method. 



14 


More extensive sample calculations will have to be studied to supple- 
ment the preliminary calculations of this paper. Of particular interest 
will be a more thorough investigation of the minimum resistance obtained 
at the high amplitude discrete tone in figure 8(b). 


CONCLUSIONS 

A model has been presented which can be used to calculate acoustic 
impedance of a perforated plate liner which is excited by multiple fre- 
quencies. The nonlinear interaction effects of this spectral excitation 
can be calculated simultaneously with the effects of grazing flow velocity 
and orifice fluid velocity and displacement. The results of the calcula- 
tions have shown the following. 

1. When the model is reduced to single frequency excitation, it pro- 
vides acceptable agreement with experimental data for variations in sound 
pressure level, grazing flow, and frequency. 

2. When the model is used with two frequency excitation a strong 
interaction of the amplitude at one frequency upon the resistance at the 
other frequency is observed. This effect is not as great as a very 
limited amount of experimental data would indicate. Final judgment can- 
not be made as to which is correct until more data are thoroughly in- 
vestigated. 


3. Sample impedance calculations were made with a strong discrete 
tone superimposed upon a simulated white noise spectrum. The acoustic 
resistances at all frequencies were found to increase dramatically with 
increasing discrete tone sound pressure level when there was no grazing 
flow. However, with the presence of a 400 feet per second grazing flow, 
increasing discrete tone amplitude had only a small effect upon acoustic 
resistance. 


4. With spectral excitation the acoustic resistance was found to be 
airl y uniform with frequency. A possible exception was observed at the 

frequency of the strong discrete tone where the resistance appears to be 
lower than at the other frequencies. 

5. The choice of phase relationships between the several frequencies 
of excitation appears to have only a small effect on the impedance cal- ~ 
culations . 
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Figure 2. - Steady flow specific resistance in the orifice 
as a function of orifice velocity (v 0 ) and grazing flow 
velocity (Voq). Data obtained by Pratt & Whitney and 
the Boeing Company. 
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Figure 3. - Specific acoustic resistance, comparison of 
theory and experiment at constant sound pressure 
level (160 dB). 
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Figure 4. - Specific acoustic resistance, comparison of 
theory and experiment at varying sound pressure level. 


Figure 5. - Specific acoustic reactance, comparison of 
theory and experiment at constant sound pressure 
level (160 d BE 
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Figure 6. - Specific acoustic resistance, comparison of 
theory and experiment with varying grazing flow 
velocity, v^. 
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Figure 7. - Comparison of theory and experiment for two 
frequency resistance measurements. 
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(B) CALCULATED EFECT OF VARYING AMPLITUDE OF A 
SINGLE FREQUENCY ON RESISTANCE WITH SPECTRAL 
EXCITATION AND 0. 
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Figure 9. - Calculated effect of varying amplitude of a 
single frequency on acoustic resistance with spectral 
excitation and * 400 ft/ sec. 
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Figure 10. - Calculated variation of acoustic resistance 
with choice of phase angles for spectral pressure 
components. 


NASA-Lewit-Com’ 




